TWO-TERM RECURRENCE FORMULAE FOR 
INDEFINITE ALGEBRAIC INTEGRALS 



DETMAR MARTIN WELZ 

Abstract. Two-term recurrence relations are supplied for indefinite integrals 
of functions that involve factors of the types P 2 ™, P 3 ™, P 4 n , Pi m Qi n , E 1 P 1 n , 
Pi m Q2 n , E 1 P 2 n , P 2 m Q 2 n , Pi m Qi n 5i", £iPi m Qi™, P 1 m Q 1 n S 2 p , and 
Pi m Qi n Si p T\ q , where Pi, Qj, and TJ denote arbitrary polynomials of 
degree i, j, k and I in the integration variable, E\ represents the exponential 
function of an arbitrary linear polynomial in this variable, and m, n, p and 
q are arbitrary constant exponents. The 136 relations leave the form of an 
integrand unchanged and increment or decrement the exponents in steps of 
unity. 



Introduction 

To be maximally useful in science and engineering applications, the symbolic 
evaluation of an indefinite integral should be expressed as compactly as possible and 
in a form suitable for verification by differentiation, and hold in the entire complex 
plane whenever feasible. With free parameters in the integrand taken to be real, 
it should further be continuous on the real line where the integrand is integrable 
along it, and preferably be real on this line if the integrand is real everywhere 
on it, and also not explicitly involve the imaginary unit unless the integrand does 
too. When these goals are in partial conflict, the resolution should be sought by 
sacrificing compactness first, absence of imaginary offset and imaginary unit next, 
and continuity at integrable poles last. 

These principles will find, if not adoption, then ready approval. Indeed, the vast 
majority of the integrands arising in practice do not involve the imaginary unit 
or complex parameters (or else they easily separate into real and imaginary parts) 
and their antiderivatives are needed for real values of the integration variable only, 
where the absence of unexpected, symbolically often hard to localize discontinuities 
ensures that definite integrals calculated according to Newton and Leibniz will be 
correct. By default, an integrator should therefore place the unavoidable branch 
cuts such in an antiderivative that they do not traverse the real line, employing 
to this end functional equations for the logarithmic and related functions impli- 
cated in those cases. Note that the details may depend on the actual function 
implementation. 

To the recurrent chagrin of countless users, the automated integrators of today's 
computer-algebra systems do not systematically respect these principles. Consider 
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in which the integrand is analytic everywhere on the real line. The first reduc- 
tion step invokes a recurrence formula to raise the exponent —2 in the integrand 
to —1, while the second step employs transformations and evaluations chosen for 
compliance with the principles stated. The major commercial contenders, Maple 
(version 16) and Mathematica (version 8), here produce discontinuous antideriva- 
tives of about ten times this size; Maple is able to correct for the discontinuity (an 
imaginary step at x = 1), whereas Mathematica bungles definite integrals that 
straddle the discontinuity (a real jump near x = —0.493). 

The short example already illustrates the potential of combining suitable re- 
currence relations with carefully chosen terminal evaluations to achieve fast and 
strong automated symbolic integration in accordance with the above principles. 
This approach is currently being explored in the rule-based integrator Rubi [3] , for 
instance. The purpose of the present publication is to furnish efficient recurrence 
formulae for the reduction of algebraic integrals that are regularly encountered in 
practice; recipes for a properly principled evaluation of the terminal instances of 
the integrals are to be published separately. 

Being two-term relations, the recurrence formulae presented here replace an al- 
gebraic integrand with a simpler one while liberating a corresponding algebraic part 
of the antiderivative. Specifically, they apply to integrands that involve products of 
arbitrary powers of linear or quadratic polynomials, or isolated arbitrary powers of 
cubic or quartic polynomials; this naturally includes integer positive and negative 
powers. The relations increment or decrement the exponents in steps of unity, ei- 
ther separately or in up-down pairs, and leave the form of an integrand unchanged. 
Accordingly, they should lend themselves especially to the automated symbolic in- 
tegration of algebraic functions in which polynomials are raised to high positive or 
negative powers. 

Also covered are integrands containing the exponential factor exp(a + bx) along 
with such powers of polynomials; this exponential derives from the power of a linear 
polynomial via Euler's famous limit: exp(a + bx) — lim„_ ! . 00 (l + a/n + b/n x) n . 
Insofar as simple analog formulae for the trigonometric factors cos(a + bx) and 
sin(a + bx) exist, they are supplied as well; otherwise these integrands should be 
split such that the integrals involve either exp[i (a + bx)] or exp[— i (a + bx)], where 
+ia and +ib replace the parameters a and b. Note that the exponential factors 
cannot be recurred on because their parameters are unrelated to the exponent of 
the linear polynomial. 

Quite generally, indefinite integrals of the present kind with generic exponents 
evaluate to special cases of those hypergeometric functions that admit one-dimen- 
sional Euler-type integral representations. With the sum of polynomial degrees 
in the algebraic integrands ranging from two to four, the Gauss function 2-F1, the 
Appell function and the Lauricella function Fp of three variables result, and 
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where an exponential factor accompanies polynomials with a combined degree of 
one or two, the incomplete gamma function (which is a special case of Kummcr's 
confluent hypergeometric function) and the Humbert function <?i appear. The 
recurrence relations for the integrals are thus equivalent to relations among con- 
tiguous instances of these hypergeometric functions, and they are of interest mainly 
for integrals whose evaluation reduces to elementary functions or to the canonical 
elliptic integrals. 

Recurrence formulae for integrands that include the power of a linear polynomial, 
say (a + bx) m , admit profitable transformations once this factor is specialized to 
the monomial x m by setting a = and b = I. Relation (|9.1.1[) from Section l9Tl 
below, for example, can thus be specialized to 

(n+l)(ad-bc) a J (A + Bx) x m (a + bx) n (c + dx) p dx 

- / [(n+l)A(ad-bc) - (m+l)(Ab - Ba) c — (m + n+p+3)(Ab — Ba) dx\ 

x rn {a + bx) n+1 (c + dx) p dx 

- (Ab - Ba) x rn+1 (a + bx) n+1 {c + dx) p+1 = 0. 

The substitution x <— x q , m<— (m+l)/q — 1 now yields 

(n+l)(ad-bc) a J (A + Bx q ) x m (a + bx q ) n (c + dx q ) p dx 

-J [(n+l)A(ad~bc) -X(Ab-Ba)c - {\ + n+p+2)(Ab - Ba) dx q ] 

x m (a + bx q ) n+1 {c + dx q ) p dx 

--(Ab- Ba) x m+l (a + bx q ) n+1 (c + dx q ) p+1 = 0, 
9 

where A = (m+l)/q. In the new kind of integrand, the monomial exponent can be 
changed only in steps of size q, and in the present relation the monomial can be 
eliminated altogether from the integrands by setting m — 0. Recurrence formulae 
of this type are needed for the reduction of elliptic integrals in particular, for which 
q = 2. 

In a further step, the monomial power may be split by substituting m 4— m + 
rn + sp, and then x rn and x sp be merged with the respective other factors: 

(n+l)(ad -bc)aj{A + Bx q ) x m {ax r + bx r+q ) n {cx s + dx s+q ) p dx 

-J [(n+1) A (ad -be)- fi {Ab - Ba)c- {ji + n+p+2){Ab - Ba) dx q ] 

x m ~ r {ax r + bx r+q ) n+1 {cx s + dx s+q ) p dx 

--(Ab- Ba) x m - r - s+1 {ax r + bx r+q ) n+1 (cx s + dx s+q ) p+1 = 0, 
9 

where p — (m + rn + sp + l)/q. The monomial exponent has now clearly lost 
its independence. Note that recurrence relations transformed in these ways hold 
for integrands involving arbitrary fractional exponents because possible branch-cut 
crossings affect all terms of a relation equally; thus m, n, p, q, r and s in the above 
are unrestricted apart from the proviso that q ^ 0. The same transformations can 
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obviously be applied to relations in which the integrand includes an exponential 
factor. The conversion of the formulae presented below, however, is left to the 
fortunate reader since this is a fairly straightforward task whose combined results 
would have consumed an unwarranted amount of space here. 

The present formula collection is kept at a comfortable size by limiting the sum 
of polynomial degrees in the algebraic integrands at four; to just include, in the 
notation of the Abstract, the degree-five product Pi" Q4" in the same way would 
roughly double the size already. Among the degree-four integrands, no combination 
is omitted but Pi" 1 Q3", as it could aid a systematic handling of elliptic integrals 
(for integer m and odd integer 2n) alongside its cousin Pi" 1 Q A n only. The combined 
degree of powers that accompany an exponential (or trigonometric) factor of linear 
argument is even limited at two: the degree-three products EiP 3 n , EiPi m Q2", and 
EiP\ m Qi" Si p already appear to be of no interest in practice, and their integrals 
do not relate to one of the established hypergeometric functions either. Exponential 
factors with higher-degree arguments are ignored for the same reason. 

Still, when adapted relations are applied to suitably transformed integrals, the 
present formulae suffice to handle any elliptic integrand. Relations for the inte- 
grand factors x m Q2" and x m Q\" S , 2 P with the substitution x <— x 2 in particular 
encompass two-term equivalents of the multiterm recurrence formulae traditionally 
employed to reduce general elliptic integrals after the x 1 and x 3 terms of the rad- 
icand have been annihilated by a Mobius transformation and the even part of the 
rational cofactor has been expanded into partial fractions with respect to x 2 . Note 
that the Mobius transformation generates an ugly piecewise-constant prefactor and 
that the expansion step may introduce complex coefficients, but these effects beset 
any reduction to canonical elliptic integrals. 

Recurrence relations for rational integrands multiplied by a power P 3 n or P4™ 
are actually made redundant for any integer An in this way; the direct reduction 
being much easier and more general, however, it is supported here for the isolated 
powers, although the terminal evaluation may still require a Mobius transforma- 
tion. Moreover, a piecewise-constant prefactor would arise for odd 2n even if an 
isolated radicand has repeated roots, but is then better avoided, as it is rarely 
a natural part of the elementary antiderivative that results. The other degree- 
three and degree-four combinations of powers covered below must be regarded as 
indispensable anyway; so the elliptic cases among them are best reduced like the 
elementary ones, and the fractional powers just collected into a single factor prior 
to the terminal evaluation, introducing a piecewise constant which fuses with that 
from a subsequent Mobius transformation. 

For each integrand, some of the recurrence relations will fail if certain polynomial 
roots coincide, because the discriminants or resultants appear as factors in front 
of the integrals. Such confluent roots can be fused, and thereby made harmless, 
through a rational refactorization of the integrand, which lowers the total degree, 
but introduces ugly piecewise-constant prefactors in general. Dedicated formulae 
are therefore provided here for all cases of degeneracy from confluent polynomial 
roots, whereby the need for refactorization is avoided, or at least deferred until 
the terminal integral is evaluated. In formulating mutually exclusive conditions for 
the various cases of degeneracy below, the polynomial factors have been assumed 
to possess nonzero leading coefficients. The conditions and associated recurrence 
formulae are most easily verified on polynomials having symbolic roots of the various 
possible confluence patterns. 
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The two-term recurrence relations have been derived by the method of undeter- 
mined coefficients and can be checked simply by differentiation. To facilitate their 
use in computer- algebra systems, they are listed here in a linear functional notation 
without subscripts or superscripts, which should be self-explaining. Extensive ab- 
breviation of repeated subexpressions on both the integrand and the formula level 
serves to condense the listing and to expose the structure of the relations. While 
two-term formulae for the simplemost integrands (a + bx + cx 2 ) n , (a + bx) m (c + dx) n 
and exp(a + bx)(c + dx) n can be found in standard tables of integrals (so in Grad- 
shteyn-Ryzhik [3] as item 2.260.2; as items 2.151, 2.153, 2.155; and for a = c = 0, 
d = 1 as items 2.321.1, 2.324.1), they are included below for convenience. 

The state of the art prior to the present publication is reflected in a 1948 book 
by Timofeyev [4] (there is no evidence for the 1933 edition cited in [2]) and a recent 
communication by Barnett pQ. 

1. Integrands involving (a + b*x + c*x~2)~n 

To exclude integrands with confluent roots, the following recurrence should be 
applied only if the quadratic discriminant does not vanish: 4*a*c-b~2^0. Note 
that the relation can be used backwards for its own inverse. 

(1) (n+ l)*(4*a*c-b~2)*INT((a+b*x + c*x~2) ~n, x) 
-2*(2*n + 3)*c*INT((a+b*x + c*x~2)~(n+l) , x) 
+ (b + 2*c*x)*(a+b*x + c*x~2)~(n+l) =0. 

2. Integrands involving (a + b*x + c*x~2 + d*x~3) ~n 
Abbreviations used with these integrands: 
ra = 2*(3*b*d- c~2), rb = 9*a*d-b*c, rc = 2* (3*a*c - b~2) . 

2.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the cubic discriminant does not 
vanish: ra*rc - rb~2 + 0. 

(2.1.1) (n+ l)*(ra*rc-rb~2)*INT((A + B*x)*(a + b*x + c*x~2 + d*x~3)~n, x) 
+ 3*INT(((3*n+4)*U*d- (2*n + 3)*V*c- (3*n + 5)*V*d*x) 

*(a + b*x + c*x~2 + d*x~3)~(n+l) , x) 

- (U*c - 2*V*b + 3* (U*d - V*c) *x - 3*V*d*x~2) 
*(a + b*x + c*x~2 + d*x~3)~(n+l) =0, 

where U = A*rb - B*rc, V = A*ra - B*rb. 

(2.1.2) 3*v*(v + d)*INT((A + B*x)*(a + b*x+ c*x~2 + d*x~3) ~n, x) 

- n*INT ( (U*rb + B*v*rc + (U*ra + B*v*rb) *x) 
*(a + b*x + c*x~2 + d*x~3)~(n- 1) , x) 

- (U*c + 2*B*v*b + 3* (U*d + B*v*c) *x + 3*B*v*d*x~2) 
*(a + b*x + c*x~2 + d*x~3)~n=0, 

where U= (3*n + 2)*A*d- (2*n+ l)*B*c, v= (3*n+ l)*d. 

2.2. Singly degenerate integrands. Here the cubic polynomial is required to 
possess a double root whence its discriminant must vanish: ra*rc - rb~2 = 0. To 
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exclude doubly degenerate cubics, the recurrence should be applied only if addi- 
tionally ra^ 0. Note that the relation can be used backwards for its own inverse. 

(2.2) (n+ l)*(2*n + 1) * (ra*c - 3*rb*d) *INT( (a+ b*x+ c*x~2 + d*x~3) ~n, x) 
+ 3*(3*n + 4)*w*d*INT((a+b*x+c*x~2 + d*x~3)~(n+l) , x) 

- (n+l)*(rb*c-2*rc*d+ (ra*c-rb*d)*x + ra*d*x~2) 
*(a + b*x + c*x~2 + d*x~3)~n 

- w* (c + 3*d*x) * (a + b*x + c*x~2 + d*x~3) ~ (n+ 1) = 0, 

where w= (3*n + 2)*d. 

3. Integrands involving (a + b*x+c*x~2 + d*x~3 + e*x~4)~n 

Abbreviations used with these integrands: 

ra = 8*c*e-3*d~2, rb = 6*b*e - c*d, rcl = 16*a*e - b*d, 

rc2 = 4*a*e + 2*b*d-c~2, rd= 6*a*d-b*c, re = 8*a*c - 3*b~2. 

sa = ra*rc2- rb~2, sb = 1/2 * (ra*rd- rb*rcl) , sc = rb*rd - rcl*rc2, 
sd = 1/2 * (rb*re - rcl*rd) , se = rc2*re - rd~2. 

3.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the quartic discriminant docs 
not vanish: sa*re - 2*sb*rd+ sc*rcl ^ 0. 

(3.1.1) (n+ 1) * (sa*re - 2*sb*rd+ sc*rcl) *INT( (A+ B*x + C*x~2) 

*(a + b*x + c*x~2 + d*x~3 + e*x~4)~n, x) 

- 4*INT ( ( (4*n+ 5) *U*e - (3*n + 4) *V*d + u*W*c 

- (2*u*V*e- (3*n + 5)*W*d)*x + (4*n + 7)*W*e*x~2) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4) " (n+ 1) , x) 

+ (U*d - 2*V*c + 3*W*b + 4* (U*e - V*d + W*c) *x 

- 4* (V*e - W*d) *x~2 + 4*W*e*x~3) 

*(a + b*x+ c*x~2 + d*x~3 + e*x~4) ~(n+ 1) =0, 

where U = A*sc - B*sd + C*se, V = A*sb - B*sc + C*sd, W = A*sa-B*sb + C*sc, 
u = 2*n + 3. 

(3.1.2) 4* (2*n + 1) * (4*n+ 3) *w*e~2*INT( (A + B*x + C*x~2) 

*(a + b*x + c*x~2 + d*x~3 + e*x~4)~n, x) 

- n*INT ( ( (U- V*v) *rcl + (V*rd + W*re) *w+ 2* ( (U - V*v) *rb 

+ (V*rc2 + W*rd)*w)*x + ((U-V*v)*ra + (V*rb + W*rcl) *w) *x~2) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4)~(n- 1) , x) 

- ((U-V*v)*d + (V*c + 3*W*b)*w + 4*((U + n*V*d)*e + W*w*c)*x 
+ 2* (V*e + 2*W*d) *w*x~2 + 4*W*w*e*x~3) 

*(a + b*x + c*x~2 + d*x~3 + e*x~4)~n=0, 

where U = 2*(2*n + 1) * ( (4*n+ 3) *A*e - (2*n + l)*C*c)*e, V= (4*n + 3)*B*e- 
(3*n+2)*C*d, W= (2*n+l)*C*e, v= (3*n+l)*d, w = 2*(4*n+ l)*e. 

3.2. Singly degenerate integrands. Here the quartic polynomial is required to 
possess one double root whence its discriminant must vanish: sa*re - 2*sb*rd + 
sc*rcl = 0. To exclude doubly degenerate quartics, the recurrences should be ap- 
plied only if additionally sa^ 0. 
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(3.2.1) u*(u+ l)*(sa*c-3*sb*d+ 6*sc*e) *sa*INT( (A+ B*x) 

*(a + b*x+c*x~2 + d*x~3 + e*x~4)~n, x) 

- 2*INT ( (2* (4*n+ 5) * (U*rb - V) *e - (3*n + 4) * (U*ra - W) *d 

- 2* (2*n+ 3) * (U*ra - W) *e*x) 

*(a + b*x+c*x~2 + d*x~3 + e*x~4)~(n+l) , x) 

- (u+ l)*(A*sa-B*sb)*(sb*c-2*sc*d + 3*sd*e 

+ (sa*c - sb*d + sc*e) *x+ (sa*d- sb*e) *x~2 + sa*e*x~3) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4) ~n 
+ ( (U*rb - V) *d - (U*ra - W) *c + 2* (2* (U*rb - V) *e - (U*ra - W) *d) *x 

- 2*(U*ra- W)*e*x~2)*(a+b*x+ c*x~2 + d*x~3 + e*x~4) ~(n+ 1) =0, 

where U= (3*n + 2)*A*sa*d- (2*n+ l)*(2*A*sb*e + B*sa*c) + 

B*((3*n+ l)*sb*d- (4*n+ l)*sc*e), V= (4*n + 3)*(2*A*rc2-B*rd)*sa*e, 

W= (4*n + 3)*(2*A*rb-B*rcl)*sa*e, u = 2*n+ 1. 

(3.2.2) 8*(2*n+ l)*(4*n+ 1) *w*e~2*INT( (A+ B*x) 

*(a + b*x + c*x~2 + d*x~3 + e*x~4)~n, x) 
+ n*INT(((U*ra + V*rb)*((2*n- l)*sa*c- (3*n-2)*sb*d 
+ (4*n-3)*sc*e) - (U*rcl +V*rd)*w+ ((U*ra + V*rb) 
*((3*n- l)*sa*d-2*(2*n- l)*sb*e) - 2* (U*rb + V*rc2) *w) *x) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4)~(n- 1) , x) 

- n* (U*ra + V*rb) * (sb*c - 2*sc*d+ 3*sd*e 

+ (sa*c - sb*d + sc*e)*x+ (sa*d- sb*e)*x~2 + sa*e*x~3) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4) ~ (n- 1) 

- w* (U*d+ V*c + 2* (2*U*e + V*d) *x + 2*V*e*x~2) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4)~n=0, 

where U = 2*(2*n + l)*A*e- (3*n+ l)*B*d, V = 2*(4*n+ l)*B*e, 
w= (4*n- l)*sa*e. 

3.3. Doubly degenerate integrands. Now the quartic polynomial is required to 
possess a triple root whence the following must vanish: 6*ra*a- 3*rb*b + rcl*c = 
rcl-rc2 = 0. To exclude triply degenerate quartics, the recurrence should be ap- 
plied only if additionally ra^ 0. Note that the relation can be used backwards for 
its own inverse. 

(3.3.1) (n + 1) * (3*n + 1) * (3*n + 2) *ra~2 

*INT ( (a + b*x + c*x~2 + d*x~3 + e*x~4) ~n , x) 
+ 4*(4*n+5)*w*e*INT((a+b*x + c*x~2 + d*x~3 + e*x~4) ~ (n+ 1) , x) 

- (n + 1) * (3* (5*n+ 3) *ra*b*e - v*c + ( (22*n + 13) *ra*c*e - 3*v*d) *x 
+ 3* ( (7*n + 4) *ra*d - 2*v) *e*x~2 + 6* (2*n + 1) *ra*e~2*x~3) 

*(a + b*x + c*x~2 + d*x~3 + e*x~4) ~n 

- w*(d + 4*e*x)*(a + b*x + c*x~2 + d*x~3 + e*x~4) ~ (n+ 1) =0, 

where v= (3*n + 2)*ra*d- (4*n + 3)*rb*e, w = 6*(2*n+ l)*(4*n + 3)*e~2. 

Now the quartic polynomial is required to possess two double roots whence 
the following must vanish: ra*d-4*rb*e = rb*d-4*rc2*e = 0. To exclude triply 
degenerate quartics, the recurrence should be applied only if additionally ra^O. 
Note that the relation can be used backwards for its own inverse. 
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(3.3.2) 2*(n+ l)*(2*n+l)*(rc2-rcl) 

*INT ( (a + b*x + c*x~2 + d*x~3 + e*x~4) ~n , x) 
+ 4*w* (w+ 2) *e*INT ( (a + b*x + c*x~2 + d*x~3 + e*x~4) ~ (n+ 1) , x) 

- (n+l)*(rd+ (rcl + 2*rc2)*x + 3*rb*x~2 + ra*x~3) 
*(a + b*x + c*x~2 + d*x~3 + e*x~4) ~n 

- w*(d + 4*e*x)*(a + b*x + c*x~2 + d*x~3 + e*x~4) ~(n+ 1) =0, 

where w = 4*n + 3. 

4. Integrands involving (a + b*x)~m*(c + d*x)~n 

To exclude integrands with confluent roots, the following recurrences should be 
applied only if the resultant of the linear polynomials does not vanish: a*d - b*c ^ 0. 
Note that each relation can be used backwards for its own inverse. 

(4.1) (m+ l)*(a*d-b*c)*INT((a+b*x)~m*(c + d*x)~n, x) 
-d*(m + n + 2)*INT((a+b*x)~(m+ 1) * (c + d*x) ~n, x) 
+ (a + b*x)~(m+l)*(c + d*x)~(n+l) =0. 

(4.2) (m+ l)*b*INT((a+b*x)~m*(c + d*x)~n, x) 

+ n*d*INT((a+b*x) ~ (m+ l)*(c + d*x) ~(n- 1) , x) 
- (a + b*x)~(m+l)*(c + d*x)~n=0. 

5. Integrands involving exp(a + b*x)*(c + d*x)~n 

Note that each recurrence can be used backwards for its own (or its partner's) 
inverse. 

(5.1) (n+ l)*d*INT(exp(a+b*x)*(c + d*x)~n, x) 
+ b*INT(exp(a + b*x)*(c + d*x)~(n+l) , x) 
-exp(a + b*x)*(c + d*x)~(n+l) =0. 

(5.2) (n+ l)*d*INT(cos(a+b*x)*(c + d*x)~n, x) 
-b*INT(sin(a+b*x)*(c + d*x)~(n+l) , x) 
-cos(a + b*x)*(c + d*x)~(n+l) =0. 

(5.3) (n+ l)*d*INT(sin(a+b*x)*(c + d*x)~n, x) 
+ b*INT(cos(a+b*x)*(c + d*x)"(n+ 1) , x) 
-sin(a + b*x)*(c + d*x)~(n+l) =0. 

6. Integrands involving (a + b*x)~m*(c + d*x + e*x~2)~n 

Abbreviations used with these integrands: 

ra=2*a*e-b*d, rb = a*d - 2*b*c. 

6.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the overall discriminant does 
not vanish: (ra*a- rb*b) * (4*c*e - d~2) ^0. 

(6.1.1) (m+ l)*(ra*a-rb*b)*INT((A + B*x) 

*(a + b*x)~m*(c + d*x + e*x~2)~n, x) 

- INT( ((m+ l)*(A*ra-B*rb) - v*V*d- 2*v*V*e*x) 
*(a + b*x)~(m+ l)*(c + d*x + e*x~2)~n, x) 

-2*V*(a + b*x)~(m+ 1) * (c + d*x + e*x~2) ~ (n+ 1) =0, 

where V = A*b - B*a, v = m + 2*n + 3. 
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(6.1.2) 2*u*e*INT((A + B*x)*(a + b*x) ~m * (c + d*x + e*x~2) ~n, x) 

- INT((u*U*a + m*B*rb + (u*U*b + m*B*ra)*x) 
*(a + b*x)~(m- l)*(c + d*x + e*x~2)~n, x) 

- 2*B* (a + b*x) ~m* (c + d*x + e*x~2) " (n+ 1) = 0, 

where U = 2*A*e - B*d, u = m + 2*n + 2. 

(6.1.3) (n+1)* l/2*(ra*a-rb*b)*(4*c*e-d~2)*INT((A + B*x) 

*(a + b*x)~m*(c + d*x + e*x~2)~n, x) 
+ INT((U*u- V*((n+ l)*ra+ a*e) - V*(u + b)*e*x) 
*(a + b*x)~m*(c + d*x + e*x~2)~(n + 1) , x) 

- (U- V*e*x)*(a + b*x) ~ (m+ 1) * (c + d*x + e*x~2) " (n+ 1) =0, 

where U = A*rb*e- (A*d - B*c) *ra, V = A*ra-B*rb, u= (m + 2*n + 3)*b. 

(6.1.4) 2*v*(v + b)*e*INT((A + B*x)*(a + b*x) ~m * (c + d*x + e*x~2) ~n, x) 
+ n*INT ( (U*rb - V*v*a + (U*ra - V*v*b) *x) 

*(a + b*x)~m*(c + d*x + e*x~2)~(n- 1) , x) 

- (U + B*v*d + 2*B*v*e*x) * (a+ b*x) " (m+ 1) * (c + d*x + e*x~2) ~n= 0, 

where U = 2*(m + 2*n + 2)*(A*b-B*a)*e + (m+ l)*B*ra, V = B* (4*c*e - d~2) , 
v= (m + 2*n + l)*b. 

(6.1.5) (n+ l)*(4*c*e-d~2)*INT((A + B*x) 

*(a + b*x)~m*(c + d*x + e*x~2)~n, x) 

- INT ( (m*U*b + v*V*a + (v + m) *V*b*x) 

*(a + b*x)~(m- l)*(c + d*x + e*x~2)~(n + 1) , x) 
+ (U + V*x)*(a+b*x)~m*(c + d*x + e*x~2)~(n+l) =0, 

where U = A*d - 2*B*c, V = 2*A*e - B*d, v = 2*n + 3. 

(6.1.6) (m+ l)*u*b~2*INT((A+B*x)*(a + b*x) ~m * (c + d*x + e*x~2) ~n, x) 
+ n*INT ( (u*U*d + V*rb + (2*u*U*e + V*ra) *x) 

*(a + b*x)~(m + l)*(c + d*x + e*x~2)~(n- 1) , x) 

- (u*U + V*a + V*b*x) * (a + b*x) ~ (m+ 1) * (c + d*x + e*x~2) ~n= 0, 

where U = A*b-B*a, V= (m+ 1)*B, u = m + 2*n + 2. 

6.2. Singly degenerate integrands. Now the root of the linear factor is required 
to coincide with one root of the quadratic whence the corresponding resultant must 
vanish: ra*a-rb*b = 0. To exclude additional degeneracies, the recurrences should 
be applied only if also ra^O. Note that each relation can be used backwards for 
its own inverse. 

(6.2.1) (m + n+ 1) *ra*INT( (a + b*x) ~m * (c + d*x + e*x~2) ~n, x) 

- (m + 2*n + 2)*e*INT((a + b*x) " (m+ l)*(c + d*x + e*x~2) ~n, x) 
+ b* (a + b*x) ~m * (c + d*x + e*x~2) ~ (n + 1) = 0. 

(6.2.2) (m + n+ l)*v*(4*c*e-d~2)*INT((a+b*x) ~m* (c + d*x + e*x~2) ~n, x) 

- w*(w+ l)*b*e*INT((a+b*x)~m*(c + d*x + e*x~2) " (n+ 1) , x) 

+ (m*a*e + v*d + w*b*e*x) * (a+ b*x) ~m * (c + d*x + e*x~2) " (n+ 1) = 0, 

where v= (n+ l)*b, w = m + 2*n + 2. 

(6.2.3) (n+1) *ra*INT ( (a + b*x) "m * (c + d*x + e*x~2) ~n , x) 

+ (m + 2*n + 2)*b~2*INT((a+b*x)~(m- 1) * (c + d*x + e*x~2) " (n+ 1) , x) 
-b*(a + b*x)~m*(c + d*x + e*x~2)~(n + 1) =0. 
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Now the quadratic polynomial is required to possess a double root whence its 
discriminant must vanish: 4*c*e-d~2 = 0. To exclude additional degeneracies, the 
recurrences should be applied only if also ra^O. Note that each relation can be 
used backwards for its own inverse. 

(6.2.4) (m + 1) *ra*INT ( (a + b*x) ~m * (c + d*x + e*x~2) ~n , x) 
-2*(m+2*n + 2)*e*INT((a + b*x) ~ (m+ l)*(c + d*x + e*x~2) ~n, x) 
+ (d + 2*e*x)*(a + b*x)~(m + 1) * (c + d*x + e*x~2) ~n = 0. 

(6.2.5) (n+ l)*(2*n+ l)*(ra*a-rb*b) 

*INT((a + b*x)~m*(c + d*x + e*x~2)~n, x) 
-w>Kw + b)*INT((a + b*x)~m*(c + d*x + e*x~2)~(n + 1) , x) 

- (n+ l)*(rb + ra*x)*(a + b*x)~(m+ 1) * (c + d*x + e*x~2) ~n 
+ w*(a + b*x)~(m + 1) * (c + d*x + e*x~2) ~ (n + 1) =0, 

where w = (m + 2*n + 2) *b. 

(6.2.6) (n+ l)*(2*n+ 1) *ra*INT( (a+ b*x) ~m*(c + d*x + e*x~2)~n, x) 
- m* (m + 2*n + 2) *b"2 

*INT((a + b*x) ~ (m- l)*(c + d*x + e*x~2) ~ (n+ 1) , x) 

- (n + 1) * (rb + ra*x) * (a + b*x) ~m * (c + d*x + e*x~2) ~n 
+ m*b* (a + b*x) ~m* (c + d*x + e*x~2) " (n+ 1) = 0. 

7. Integrands involving exp(a + b*x)*(c + d*x + e*x~2)~n 

7.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the quadratic discriminant does 
not vanish: 4*c*e-d~2^0. 

(7.1.1) (n+ l)*(4*c*e-d~2)*INT((A + B*x) 

*exp(a + b*x)*(c + d*x + e*x~2)~n, x) 

- INT((U*b+ (2*n + 3)*V + V*b*x) 

*exp(a + b*x)*(c + d*x + e*x~2)~(n+l) , x) 
+ (U + V*x)*exp(a + b*x)*(c + d*x + e*x~2)~(n + 1) =0, 

where U = A*d - 2*B*c, V = 2*A*e - B*d. 

(7.1.2) b~2*INT((A + B*x)*exp(a+b*x)*(c + d*x + e*x~2)~n, x) 
+ n*INT( (U*d- 2*B*b*c + (2*U*e - B*b*d) *x) 

*exp(a + b*x)*(c + d*x + e*x~2)~(n- 1) , x) 

- (U + B*b*x) *exp(a+ b*x) * (c + d*x + e*x~2) ~n= 0, 

where U = A*b- (2*n+ 1)*B. 

7.2. Degenerate integrands. Now the quadratic polynomial is required to pos- 
sess a double root whence its discriminant must vanish: 4*c*e - d~2 = 0. Note that 
each recurrence can be used backwards for its own inverse. 

(7.2.1) 2*(n + l)*(2*n+ 1) *e*INT(exp(a+b*x) * (c + d*x + e*x~2) ~n, x) 
-b~2*INT(exp(a+b*x)*(c + d*x + e*x~2)~(n+l) , x) 

- (n+ l)*(d + 2*e*x)*exp(a+b*x)*(c + d*x + e*x~2)~n 
+ b*exp(a + b*x)*(c + d*x + e*x~2)~(n+ 1) =0. 
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(7.2.2) 2*(n+ l)*(2*n+ 1) *e*INT(cos (a+b*x) * (c + d*x + e*x~2) ~n, x) 
+ b~2*INT(cos(a+b*x)*(c + d*x + e*x~2)~(n+l) , x) 

- (n+ l)*(d+2*e*x)*cos(a + b*x)*(c + d*x + e*x~2)~n 
-b*sin(a + b*x)*(c + d*x + e*x~2)~(n + 1) =0. 

(7.2.3) 2*(n+ l)*(2*n + l)*e*INT(sin(a+b*x)*(c + d*x + e*x~2) "n, x) 
+ b~2*INT(sin(a+b*x)*(c + d*x + e*x~2)~(n+l) , x) 

- (n+ l)*(d+2*e*x)*sin(a + b*x)*(c + d*x + e*x~2)~n 
+ b*cos(a + b*x)*(c + d*x + e*x~2)~(n + 1) =0. 

8. Integrands involving (a + b*x+c*x~2)~m*(d + e*x + f*x~2)~n 
Abbreviations used with these integrands: 

ra = b*f - c*e, rb = a*f - c*d, rc = a*e - b*d. 

sa = ra*b - 2*rb*c, sb = ra*a- rc*c, sc = 2*rb*a- rc*b, 
sd = 2*rb*f - ra*e, se = rc*f - ra*d. 

8.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the overall discriminant docs 
not vanish: (ra*rc - rb~2) * (4*a*c - b~2) * (4*d*f - e~2) # 0. 

(8.1.1) u*(4*a*c-b~2)*INT((A+B*x + C*x~2) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~n, x) 

- INT(((2*A*c-B*b + 2*C*a)*u-v*V*e + W*((m + l)*rb- c*d) 

- (2*v*V*f - W*((m+ l)*ra- (v+ l)*c*e))*x- (2*v + 1) *W*c*f *x~2) 
*(a + b*x + c*x~2)~ (m + l)*(d + e*x + f *x~2) "n, x) 

- (V + W*c*x)*(a + b*x+c*x~2)~(m + l)*(d + e*x + f *x~2)~(n+ 1) =0, 

where V= (A*b-B*a)*sa- (A*c - C*a) *sb, W = A*sa-B*sb + C*sc, 
u= (m+ l)*(ra*rc- rb~2) , v = m + n + 2. 

(8.1.2) 2*w*(2*w+f )*c*INT((A+B*x + C*x~2) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~n, x) 

- INT ( (U*a + m*V*rc + W*w*d + (U*b + 2*m*V*rb + W*w*e) *x 

+ (U*c+m*V*ra + W*w*f )*x~2) 

*(a + b*x + c*x~2)~(m- 1) * (d + e*x + f *x~2) ~n, x) 

- (V + C*w*b + 2*C*w*c*x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~(n + 1) =0, 

where U= (m + n+ l)*((2*m+2*n + 3)*(2*A*f - B*e + 2*C*d) *f - 
(2*m+n + 2)*C*(4*d*f -e~2))*c, V= (2*m + 2*n + 3) * (B*f - C*e) *c - 
(n+ l)*C*ra, W = m*C*(4*a*c-b~2), w = (m + n+ l)*f . 

(8.1.3) u*c*INT((A + B*x + C*x~2)*(a + b*x + c*x~2) ~m*(d + e*x + f *x~2) "n, x) 

- INT ( (n*V*e + (C*u + w*W) *d + (2*n*V*f + (C*u + (w + n) *W) *e) *x 

+ (C*u + (w + 2*n)*W)*f*x~2) 

*(a + b*x + c*x~2)~ (m+ l)*(d + e*x + f *x~2) ~(n- 1) , x) 
+ (V + W*x)*(a + b*x + c*x~2)~(m+ 1) * (d+ e*x + f *x~2) ~n = 0, 

where V = A*b*c - 2*B*a*c + C*a*b, W = 2*(A*c-C*a)*c- (B*c-C*b)*b, 
u= (m+ l)*(4*a*c-b~2), w = 2*m + 3. 
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8.2. Integrands with a single cross-degeneracy. Here the two quadratics are 
required to have one root in common whence their resultant must vanish: ra*rc - 
rb~2 = 0. To exclude additional degeneracies, the recurrences should be applied 
only if also (2*a*f - b*e + 2*c*d) *ra^ 0. 

(8.2.1) (m + n+ 1) *v*sd* (4*a*c - b~2) *INT( (A + B*x) 

*(a + b*x+ c*x~2) ~m*(d + e*x + f *x~2) ~n, x) 
+ INT( (V* (u*c*e - v) + W* (u*b*f + n*ra) + 2*u* (V + W) *c*f *x) 

*(a + b*x+c*x~2)~(m+ 1) * (d + e*x + f *x~2) ~n, x) 
+ ra*(n*A*se*c+ ((m+ l)*A*b- (2*m + n + 2)*B*a)*sd + V*x) 

*(a + b*x + c*x~2) " (m+ l)*(d + e*x + f *x~2) ~n 

- (V + W)*f *(a+b*x + c*x~2) ~ (m+ 2) * (d+ e*x + f *x~2) ~n=0, 

where V= ((2*m + n + 2)*A*c- (m + n+ l)*B*b)*sd + n*B*se*c, 
W= (m+ l)*(A*sa-B*sb)*f , u = m + n + 2, v = (m+ l)*ra. 

(8.2.2) w*(w + f )*ra*c*INT((A+B*x) 

*(a + b*x+ c*x~2) ~m*(d+ e*x + f *x~2)~n, x) 

- INT((U*((m + n)*sd*b-n*se*c) +m*V*sb 
+ ((2*m + n)*U*sd*c + m*V*sa)*x) 

*(a + b*x+c*x~2)~(m- 1) * (d + e*x + f *x~2) ~n, x) 
+ (U*sd - V*rb - V*ra*x) * (a + b*x + c*x~2) ~m * (d + e*x + f *x~2) ~n 

- B*w*ra*f * (a + b*x + c*x~2) " (m+ 1) * (d + e*x + f *x~2) ~n= 0, 

where U= (m + n+ 1) * (2*A*f - B*e) *c + m*B*ra, 

V = ((m + n+ l)*(2*A*c-B*b)*f -n*B*ra)*f , w= (2*m + 2*n + l)*f . 

(8.2.3) 2*(m+ l)*sa*c*INT((A+B*x) 

*(a + b*x + c*x~2) ~m*(d+ e*x + f *x~2) ~n, x) 
+ INT( ( (m+ 1) *U*sd- V*e - 2*V*f *x) 

*(a + b*x+c*x~2)~(m+ 1) * (d + e*x + f *x~2) " (n- 1) , x) 
+ (U*ra-B*sa)*(a + b*x+ c*x~2) ~ (m+ l)*(d + e*x + f *x~2) ~n= 0, 

where U = 2*A*c-B*b, V= (m + n+ l)*(2*A*c-B*b)*ra-n*B*sa. 

8.3. Integrands with a single self-degeneracy. Here the first quadratic poly- 
nomial is required to possess a double root whence its discriminant must vanish: 
4*a*c-b~2 = 0. To exclude additional degeneracies, the recurrences should be ap- 
plied only if also (2*a*f - b*e + 2*c*d) * (4*d*f - e~2) + 0. 

(8.3.1) (m+l)*u*v-2*INT((A+B*x) 

*(a + b*x+ c*x~2) ~m*(d+ e*x + f *x~2) ~n, x) 
+ INT(((u+ l)*((n + l)*A*v-2*(u + n + 1) * (A*rb - B*rc) ) *f 
+ (u + n + 2)*W*e + 2*(m + n + 2)*W*f*x) 
*(a + b*x + c*x~2) " (m+ 1) * (d + e*x + f *x~2) ~n, x) 

- (m+ l)*v*(A*b- 2*B*a+ (2*A*c - B*b) *x) 
*(a + b*x + c*x~2) ~m*(d+ e*x + f *x~2) ~ (n+ 1) 

- W*(a + b*x + c*x~2)~(m+ 1) * (d+ e*x + f *x~2) ~ (n+ 1) =0, 

where W = 2*(2*m + n + 2)*(A*ra-B*rb) - (n+ 1) *B* (2*a*f - b*e + 2*c*d) , 
u = 2*m+ 1, v = 2*a*f - b*e + 2*c*d. 
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(8.3.2) 2*w*(w+f )*INT((A + B*x)*(a + b*x+c*x~2)~m*(d + e*x + f*x~2)~n, x) 

- INT ( (U* (u*rb + v) + u*V*r c + (u*U*ra + V* (u*rb - v) ) *x) 
*(a + b*x + c*x~2)~(m- 1) * (d + e*x + f *x~2) ~n, x) 

- (U*b + 2*V*a+ (2*U*c+V*b)*x) 

*(a + b*x + c*x~2)~ (m- l)*(d + e*x + f *x~2) ~(n+ 1) 

- 2*B*w* (a + b*x + c*x~2) ~m * (d + e*x + f *x~2) ~ (n+ 1) = 0, 

where U= (m + n+ 1) * (2*A*f - B*e) -m*B*e, V = 2*m*B*f, u = 2* (2*m + n) , 
v= (n+ l)*(2*a*f -b*e + 2*c*d), w= (2*m+2*n + l)*f . 

(8.3.3) (n+ l)*u*v*INT((A + B*x)*(a + b*x + c*x~2) ~m*(d + e*x + f *x~2) ~n, x) 

- INT ( (m* (2*A*c - B*b) *u + (2*n + 3) *V*v - w*W*e - 2*w*W*f *x) 
*(a + b*x + c*x~2)~m*(d + e*x + f*x~2)~(n+l) , x) 

+ v*(A*e- 2*B*d+ V*x)*(a + b*x + c*x~2) ~m*(d + e*x + f *x~2) ~ (n+ 1) 

- W*(a + b*x+ c*x~2) ~m*(d + e*x + f *x~2) ~ (n+2) =0, 

where V = 2*A*f -B*e, W = 2* (A*e - 2*B*d) *c - (2*A*f - B*e) *b, u = 4*d*f -e~2, 
v = 2*a*f -b*e + 2*c*d, w = m + n + 2. 

(8.3.4) 2*w*(w+ l)*c*f *INT((A+B*x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~n, x) 
+ n*INT((U*e- 2*V*d + (2*U*f - V*e) *x) 

*(a + b*x + c*x~2)~m*(d + e*x + f*x~2)~(n- 1) , x) 

- (U + V*x) * (a + b*x + c*x~2) ~m * (d + e*x + f *x~2) "n 

- 2*B*w*c* (a + b*x + c*x~2) ~m*(d + e*x + f *x~2) ~(n+ 1) =0, 

where U= (m + n+ l)*(2*A*f -B*e)*b + B*(2*m*rb- (m + 2*n+l)* 
(2*a*f - b*e + 2*c*d) ),V = 2*((m + n+ l)* (2*A*f - B*e) *c + m*B*ra) , 
w = 2*m + 2*n+ 1. 

(8.3.5) 2*(m+ l)*(2*m + l)*(2*a*f -b*e + 2*c*d)*c*INT((A+B*x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~n, x) 
+ n*INT( (W* (4*d*f - e~2) + V*e + 2*V*f *x) 

*(a + b*x + c*x~2)~ (m + l)*(d + e*x + f *x~2) ~(n- 1) , x) 

- (V-W*e-2*W*f *x)*(a+b*x + c*x~2)~(m+ 1) * (d+ e*x + f *x~2) ~n 
-W*(b + 2*c*x)*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~(n+ 1) =0, 

where V = 2*(m + n + 1) * (A*ra- B*rb) + (m - n) *B* (2*a*f - b*e + 2*c*d) , 
W= (m+l)*(2*A*c-B*b). 

(8.3.6) 2*(n + l)*u*f*INT((A+B*x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~n, x) 
+ INT ( (m*B*u*b - V* (v*rb - (m - n - 2) * (2*a*f - b*e + 2*c*d) ) 
+ (2*m*B*u*c-v*V*ra)*x) 

*(a + b*x + c*x~2)~ (m- l)*(d + e*x + f *x~2) ~(n+ 1) , x) 
+ 2*f *(A*e- 2*B*d + V*x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~(n + 1) 
-V*(b + 2*c*x)*(a + b*x + c*x~2)~(m- 1) * (d + e*x + f *x~2) ~ (n+ 2) =0, 

where V = 2*A*f -B*e, u = 4*d*f -e~2, v = 2*(m + n + 1) . 
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8.4. Doubly degenerate integrands. Now the two quadratics are required to 
have two roots in common whence the following must vanish: ra=rb = 0. To 
exclude additional degeneracies, the recurrences should be applied only if also 
2*a*f - b*e + 2*c*d^ 0. Note that each relation can be used backwards for its 
own inverse. 

(8.4.1) u*(4*a*c-b~2)*INT((a+b*x+ c*x~2) ~m*(d + e*x + f *x~2) ~n, x) 

- 2*(2*u + l)*c*INT((a + b*x + c*x~2) ~ (m+ l)*(d + e*x + f *x~2)~n, x) 
+ (b + 2*c*x)*(a + b*x + c*x~2)~(m + l)*(d + e*x + f *x~2)~n=0, 

where u = m + n + 1 . 

(8.4.2) c*INT((a + b*x+c*x~2)~m*(d + e*x + f*x~2)~n, x) 

-f *INT((a + b*x + c*x~2)~ (m+ 1) * (d+ e*x + f *x~2) ~(n- 1) , x) =0. 

Now the first quadratic is required to possess a double root which coincides with 
one root of the second quadratic whence the following must vanish: 4*a*c-b~2 = 
2*a*f - b*e + 2*c*d = 0. To exclude additional degeneracies, the recurrences should 
be applied only if also 4*d*f -e~2 + 0. Note that each relation can be used back- 
wards for its own inverse. 

(8.4.3) (2*m + n+ l)*v*(4*d*f - e~2)*c 

*INT((a + b*x + c*x~2) ~m*(d + e*x + f *x~2) "n, x) 
+ 2*w*(2*w+f )*INT((a + b*x + c*x~2) " (m+ l)*(d + e*x + f *x~2)~n, x) 

- (v*ra + w*b + 2*w*c*x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~(n+ 1) =0, 

where v = 2*m + n + 2, w = (m + n+ l)*f . 

(8.4.4) (2*m + n + l)*(n + l)*(4*d*f - e~2)*c 

*INT((a + b*x + c*x~2) ~m*(d + e*x + f *x~2) ~n, x) 
-2*w*(2*w+l)*c*f 

*INT((a + b*x + c*x~2)~m*(d + e*x + f*x~2)~(n+l) , x) 
+ (m*ra + w*c*e + 2*w*c*f *x) 

*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~(n + 1) =0, 

where w = m + n + 1 . 

(8.4.5) 2*(2*m + n + 1) *c*INT( (a + b*x + c*x~2) ~m* (d + e*x + f *x~2)~n, x) 

+ 2*n*f *INT((a + b*x + c*x~2)~(m+ l)*(d+e*x + f *x~2) ~(n- 1) , x) 

- (b + 2*c*x)*(a + b*x + c*x~2) ~m*(d + e*x + f *x~2) ~n= 0. 

Now both quadratics are required to possess a double root whence the following 
must vanish: 4*a*c - b~2 = 4*d*f - e~2 = 0. To exclude additional degeneracies, the 
recurrences should be applied only if also ra^O. Note that each relation can be 
used backwards for its own inverse. 

(8.4.6) 2*(2*m+ 1) *w*c*INT( (a+ b*x + c*x~2) ~m*(d + e*x + f *x~2) ~n, x) 
-4*v*(2*v+ l)*c*f 

*INT((a + b*x + c*x~2) ~ (m+ l)*(d + e*x + f *x~2)~n, x) 

- 2* (n*ra - v*b*f - 2*v*c*f *x) 

*(a + b*x + c*x~2)~(m+ 1) * (d + e*x + f *x~2) ~n 
-w*(b + 2*c*x)*(a + b*x + c*x~2)~m*(d + e*x + f *x~2)~n = 0, 

where v = m + n+l, w = (m+ 1) * (2*a*f - b*e + 2*c*d) . 



RECURRENCE FORMULAE FOR ALGEBRAIC INTEGRALS 



15 



(8.4.7) 2*(m+ l)*(2*m+ 1) *c*INT( (a+ b*x + c*x~2) ~m*(d + e*x + f *x~2) ~n, x) 
-2*n*(2*n- l)*f 
*INT((a + b*x + c*x~2)~(m + 1) * (d+ e*x + f *x~2) " (n- 1) , x) 

- (m+ l)*(b + 2*c*x)*(a + b*x + c*x~2) ~m*(d + e*x + f *x~2) ~n 

+ n*(e + 2*f*x)*(a + b*x+c*x~2)~(m + 1) * (d + e*x + f *x~2) " (n- 1) =0. 

9. Integrands involving (a + b*x)~m*(c + d*x)~n*(e + f*x)~p 
Abbreviations used with these integrands: 
ra= (c*f + d*e)*b- a*d*f , rb = b*c*e. 

9.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the overall discriminant does 
not vanish: (a*d - b*c) * (a*f - b*e) * (c*f - d*e) i= 0. 

(9.1.1) (m+ l)=Ka*d-b*c)=Ka*f -b*e)*INT((A + B*x) 

* (a + b*x) ~m* (c + d*x) ~n* (e + f *x) ~p, x) 

+ INT(((m+ l)*(A*ra-B*rb) + V*((n + l)*d*e+ (p+ l)*c*f ) 
+ (m + n + p + 3)*V*d*f *x) 

*(a + b*x)~(m+ l)*(c + d*x)~n*(e + f*x)~p, x) 

- V* (a + b*x) " (m+ 1) * (c + d*x) " (n+ 1) * (e + f *x) " (p+ 1) = 0, 

where V = A*b - B*a. 

(9.1.2) (m + n + p + 2)*d*f*INT((A+B*x) 

* (a + b*x) ~m* (c + d*x) ~n * (e + f *x) ~p, x) 

- INT ( (U*a - m*B*rb + (U*b - m*B*ra) *x) 

*(a + b*x)~(m- l)*(c + d*x)~n*(e + f *x)~p, x) 
-B*(a + b*x)~m*(c + d*x)~(n + 1) * (e + f *x) " (p+ 1) = 0, 

where U= (m + n + p + 2)*A*d*f -B*((n+ l)*d*e + (p+ l)*c*f ). 

(9.1.3) u*b*INT ( (A + B*x) * (a+ b*x) ~m * (c + d*x) ~n * (e + f *x) ~p , x) 
- INT((B*u*c+V*(v*c*f + n*d*e) + (B*u + (v + n)*V*f )*d*x) 

*(a + b*x)~(m + l)*(c + d*x)~(n- 1) * (e + f *x) ~p, x) 
+ V*(a + b*x) ~ (m+ l)*(c + d*x) ~n * (e + f *x) " (p+ 1) = 0, 

where V = A*b-B*a, u= (m+ l)*(a*f -b*e), v = m + p + 2. 

9.2. Singly degenerate integrands. Now the roots of the first two linear fac- 
tors are required to coincide whence the corresponding resultant must vanish: 
a*d-b*c = 0. To exclude additional degeneracies, the recurrences should be ap- 
plied only if also a*f - b*e + 0. Note that each relation can be used backwards for 
its own inverse. 

(9.2.1) (m + n+ l)*(a*f -b*e)*INT((a+b*x) ~m*(c + d*x) ~n*(e + f *x) ~p, x) 

- (m + n + p + 2)*f*INT((a+b*x)~(m+ 1) * (c + d*x) ~n * (e + f *x) ~p, x) 
+ (a + b*x) ~(m + l)*(c + d*x) ~n * (e + f *x) " (p+ 1) =0. 

(9.2.2) (p+ l)*(a*f -b*e)*INT((a+b*x) ~m*(c + d*x) ~n * (e + f *x) ~p , x) 

+ (m + n + p + 2)*b*INT((a+b*x)~m*(c + d*x)~n*(e + f*x)~(p+l) , x) 

- (a + b*x) ~ (m+ l)*(c + d*x) ~n * (e + f *x) " (p+ 1) =0. 

(9.2.3) b*INT((a + b*x)~m*(c + d*x)~n*(e + f*x)~p, x) 
-d*INT((a + b*x) ~ (m+ l)*(c + d*x)~(n- l)*(e + f *x) ~p, x) =0. 
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(9.2.4) (m + n+ 1) *b*INT( (a + b*x) ~m* (c + d*x) ~n * (e + f *x) ~p, x) 

+p*f *INT((a+b*x) ~ (m+ 1) * (c + d*x) "n * (e + f *x) ~ (p- 1) , x) 
- (a + b*x) ~ (m+ l)*(c + d*x)~n*(e + f*x)~p = 0. 

10. Integrands involving exp(a + b*x)*(c + d*x)~m*(e + f*x)~n 

10.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the resultant of the linear poly- 
nomials does not vanish: c*f - d*e + 0. 

(10.1.1) u*(c*f -d*e)*INT((A + B*x) 

*exp (a + b*x) * (c + d*x) ~m * (e + f *x) ~n, x) 

- INT(((A*f -B*e)*u + V*(b*e + (n+ l)*f ) + V*b*f *x) 
*exp(a + b*x)*(c + d*x)~(m + l)*(e + f*x)~n, x) 

+ V*exp(a + b*x)*(c + d*x)~(m+ 1) * (e + f *x) " (n+ 1) =0, 

where V = A*d-B*c, u= (m+ l)*d. 

(10.1.2) b*f *INT((A + B*x)*exp(a+b*x)*(c + d*x) ~m* (e + f *x) ~n, x) 

- INT((U*c-m*B*d*e+ (U-m*B*f )*d*x) 
*exp(a + b*x)*(c + d*x)~(m- l)*(e + f*x)~n, x) 

-B*exp(a + b*x)*(c + d*x)~m*(e + f*x)~(n+ 1) =0, 

where U= (A*f -B*e)*b- (n+ l)*B*f . 

(10.1.3) u*d*INT ( (A + B*x) *exp (a+ b*x) * (c + d*x) ~m * (e + f *x) ~n, x) 

- INT((B*u*e-V*(b*e + n*f ) + (B*u - V*b) *f *x) 
*exp(a + b*x)*(c + d*x)~(m+ 1) * (e + f *x) " (n - 1) , x) 

- V*exp(a + b*x)*(c + d*x)~(m+ 1) * (e + f *x) ~n= 0, 

where V = A*d-B*c, u= (m+ l)*d. 

10.2. Degenerate integrands. Now the roots of the two linear factors arc re- 
quired to coincide whence the corresponding resultant must vanish: c*f - d*e = 0. 
Note that each recurrence can be used backwards for its own (or its partner's) 
inverse. 

(10.2.1) (m + n+ l)*d*INT(exp(a+b*x)*(c + d*x) ~m * (e + f *x) ~n, x) 
+ b*INT(exp(a+b*x)*(c + d*x) ~ (m+ l)*(e + f *x)~n, x) 
-exp(a + b*x)*(c + d*x)~(m + l)*(e + f *x)~n = 0. 

(10.2.2) (m + n+ l)*d*INT(cos(a+b*x)*(c + d*x) ~m * (e + f *x) ~n, x) 
-b*INT(sin(a+b*x)*(c + d*x) ~ (m+ l)*(e + f *x) ~n, x) 

- cos(a + b*x)*(c + d*x) ~(m + l)*(e + f *x) ~n = 0. 

(10.2.3) (m + n+ l)*d*INT(sin(a+b*x)*(c + d*x) ~m * (e + f *x) ~n, x) 
+ b*INT(cos(a+b*x)*(c + d*x) " (m+ l)*(e + f *x) ~n, x) 

- sin(a + b*x)*(c + d*x) ~ (m+ l)*(e + f *x) ~n=0. 

(10.2.4) d*INT(exp(a+b*x)*(c + d*x)~m*(e + f*x)~n, x) 

- f *INT(exp(a+b*x)*(c + d*x)~(m+ l)*(e + f *x) ~ (n- 1) , x) =0. 

(10.2.5) d*INT(cos(a+b*x)*(c + d*x)~m*(e + f*x)~n, x) 

- f *INT(cos(a+b*x)*(c + d*x)~ (m+ l)*(e + f *x) ~ (n- 1) , x) =0. 

(10.2.6) d*INT(sin(a+b*x)*(c + d*x)~m*(e + f*x)~n, x) 

- f *INT(sin(a+b*x)*(c + d*x)~(m+ l)*(e + f *x) ~(n- 1) , x) =0. 
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11. Integrands involving (a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p 
Abbreviations used with these integrands: 

ra = 2*a*g-b*f , rb = a*f - 2*b*e, rc = 2*c*g - d*f , rd = c*f - 2*d*e. 
re=l/2*(ra*d + rc*b), rf = 1/2 * (ra*c + rd*b) , rg= 1/2 * (rb*c + rd*a) . 
se = 2*rf*g-re*f, sf = rg*g- re*e, sg = rg*f - 2*rf *e. 

11.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the overall discriminant does not 
vanish: (a*d - b*c) * (ra*a- rb*b) * (rc*c - rd*d) * (4*e*g- f ~2) #0. 

(11.1.1) u*l/2*(ra*a-rb*b)*INT((A + B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ INT( ( ( (A*f - B*e) *b - U*a) *u - V* (v*e + (p + 1) *rd) 

+ ((U*b- (B*g-C*f )*a)*u- V*(v*f + (p + l)*rc))*x-V*v*g*x~2) 
*(a + b*x)~(m+ l)*(c + d*x) ~n*(e + f *x + g*x~2)~p, x) 
+ V*(a + b*x) ~ (m+ 1) * (c + d*x) " (n+ 1) * (e + f *x + g*x~2) " (p+ 1) =0, 

where U = A*g-C*e, V = A*b~2 - B*a*b + C*a~2, u= (m+ l)*(a*d-b*c), 
v= (m + n + 2*p + 4)*d. 

(11.1.2) t*g*INT((A + B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT((t*T*a + V*e- W*rg + (t*(T*b + U*a) + V*f - 2*W*rf ) *x 

+ (t*U*b + V*g - W*re) *x~2) 

*(a + b*x)~(m- l)*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- C*(a + b*x) ~m*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~ (p+ 1) =0, 

where T = A*g-C*e, U = B*g-C*f, V= (m + p + l)*C*(a*d-b*c), W= (p+l)*C, 
t= (m + n + 2*p + 3)*d. 

(11.1.3) u*(4*e*g-f ~2)*INT((A + B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT( ( (2*A*g- B*f + 2*C*e) *u + V* (v*b*c + w*a*d) 

+ W*(a*c*g+ (p+ l)*rf ) + ((v + w)*V*b*d 

+ W*(((v + l)*b*c+ (w+ l)*a*d)*g+ (p+ l)*re))*x 

+ (v + w+ l)*W*b*d*g*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f *x + g*x~2) ~ (p+ 1) , x) 
+ (V + W*g*x) 

* (a + b*x) ~ (m + 1) * (c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) " (p + 1) = 0, 

where V= (A*f - B*e) *se - (A*g- C*e)*sf , W = A*se - B*sf + C*sg, 

u= (p+ 1)* 1/2 *(ra*a-rb*b)* 1/2 * (rc*c - rd*d) , v = m + p + 2, w = n + p + 2. 

(11.1.4) 2*w* (w + b*d) *g*INT( (A+ B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ INT((V*(v*e + p*rg) - (U*e + W*a*c)*w 

+ (V*(v*f + 2*p*rf) - (U*f +W*(a*d + b*c))*w)*x 

+ (V*(v*g+p*re) - (U*g + W*b*d)*w)*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 

- (V + C*w*f + 2*C*w*g*x) 

* (a + b*x) ~ (m + 1) * (c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~p = 0, 
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where U = 2* (m + n + 2*p + 3) * (A*b*d - C*a*c) *g + C* ( (m + 1) *ra*c + 
(n + 1) *rc*a) , V = 2* (m + n + 2*p + 3) * (B*b*d - C* (a*d + b*c) ) *g + 
C*((m+ l)*ra*d + (n+l)*rc*b), W = p*C* (4*e*g- f ~2) , 
v= (m + p + l)*b*c + (n + p + l)*a*d, w= (m + n + 2*p + 2)*b*d. 

(11.1.5) u*b*INT((A + B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT((U*c- V*(v*e + w*rd) + (U*d + C*u*c - V* (v*f + w*rc) ) *x 

+ (C*u*d- V*v*g)*x~2) 

*(a + b*x)~(m+ l)*(c + d*x) ~ (n- l)*(e + f *x + g*x~2)~p, x) 

- 2*V* (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~ (p + 1) = 0, 

where U= (m+ l)*(A*ra*b- (B*b - C*a) *rb) , V = A*b~2- (B*b-C*a)*a, 
u = (m+ l)*(ra*a-rb*b), v = 2*(m + n + 2*p + 3)*d, w = m + 2*p + 3. 

(11.1.6) u*l/2 * (rc*c - rd*d) *INT( (A + B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT( (U*u*a+ V* (m*b*c + (n + 1) *a*d) + w*W*a*c + ( (U*u + v*V*d) *b 

+ W*((w + m)*b*c + (w + n+ l)*a*d))*x + (v + w)*W*b*d*x~2) 
*(a + b*x)~(m- l)*(c + d*x)~n*(e + f*x + g*x~2)~(p+ 1) , x) 
+ (V + W*x) 

*(a + b*x)~m*(c + d*x)~(n+ l)*(e + f *x + g*x~2)~(p + 1) =0, 

where U = A*d~2-B*c*d + C*c~2, V= (A*f - B*e) *rc - (A*g- C*e)*rd, 

W= (A*g-C*e)*rc- (B*g- Of )*rd, u= (p + 1) * (4*e*g- f ~2) , v = m + n+l, 

w = 2*p + 3. 

(11.1.7) v*w*b~2*INT((A + B*x + C*x~2) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT( (U*e + V*v*(u*e + p*rd) - p*C*w*rg+ (U*f + V*v*(u*f + p*rc) 

- 2*p*C*w*rf ) *x+ ( (U + V*u*v) *g - p*C*w*re) *x~2) 
*(a + b*x)~(m+l)*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 
+ (V*v - C*w*a - C*w*b*x) 

* (a + b*x) " (m+ l)*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~p = 0, 

where U= (m+ l)*((m + n + 2*p + 3)*(A*d~2-B*c*d + C*c~2)*b~2- 

(n + p + l)*C*(a*d-b*c)~2), V = A*b~2 - B*a*b + C*a~2, u = (n+ 2*p + l)*d, 

v= (m + n + 2*p + 3)*d, w= (m+ l)*(a*d-b*c). 

11.2. Integrands with linear-linear cross-degeneracy. Here the roots of the 
two linear factors are required to coincide whence the corresponding resultant must 
vanish: a*d-b*c = 0. To exclude additional degeneracies, the recurrences should 
be applied only if also (ra*a- rb*b) * (4*e*g- f ~2) # 0. 

(11.2.1) u*(ra*a-rb*b)*INT((A + B*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT( (u* (A*ra- B*rb) - v*V*f - 2*v*V*g*x) 

*(a + b*x)~(m+ l)*(c + d*x) ~n*(e + f *x + g*x~2)~p, x) 

- 2*V* (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) " (p + 1) = 0, 

where V = A*b - B*a, u = m + n+ 1, v = m + n + 2*p + 3. 
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(11.2.2) 2*u*g*INT((A+B*x) 

* (a + b*x) "m * (c + d*x) "n * (e + f *x + g*x~2) ~p , x) 

- INT((u*U*a+V*rb+ (u*U*b + V*ra) *x) 

*(a + b*x)~(m- l)*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- 2*B* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) " (p+ 1) = 0, 

where U = 2*A*g-B*f , V= (m + n)*B, u = m + n + 2*p + 2. 

(11.2.3) (p + 1) * 1/2 * (ra*a - rb*b) * (4*e*g- f ~2) *INT ( (A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
+ INT ( (U*u - V* ( (p + 1) *ra + a*g) - V* (u + b) *g*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+l) , x) 

- (U- V*g*x) 

*(a + b*x) ~ (m+ l)*(c + d*x) ~n * (e + f *x + g*x~2) " (p+ 1) =0, 
where U = A*rb*g - (A*f - B*e) *ra, V = A*ra - B*rb, u = (m + n + 2*p + 3) *b. 

(11.2.4) 2*v*(v + b)*g*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
+ p*INT ( (U*rb - V*v*a + (U*ra - V*v*b) *x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 

- (U + B*v*f + 2*B*v*g*x) 

* (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~p = 0, 

where U = 2*(m + n + 2*p + 2)*(A*b-B*a)*g + (m + n+ l)*B*ra, 
V = B*(4*e*g-f ~2), v= (m + n + 2*p + l)*b. 

(11.2.5) b*INT((A+B*x)*(a+b*x) ~m*(c + d*x) ~n * (e + f *x + g*x~2) ~p, x) 
-d*INT((A + B*x) 

*(a + b*x) ~ (m+ l)*(c + d*x) ~(n- l)*(e + f *x + g*x~2)~p, x) =0. 

(11.2.6) (p+l)*(4*e*g-f~2)*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT(((m + n)*U*b + v*V*a+ (v + m + n)*V*b*x) 

*(a + b*x)~(m- l)*(c + d*x)~n*(e + f*x + g*x~2)~(p+ 1) , x) 
+ (U + V*x)*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p + 1) =0, 

where U = A*f -2*B*e, V = 2*A*g-B*f, v = 2*p + 3. 

(11.2.7) u*v*b-2*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
+ p*INT ( (u*U*f + v*B*rb + (2*u*U*g + v*B*ra) *x) 

*(a + b*x)~(m+l)*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 

- (u*U + v*B*a + v*B*b*x) 

* (a + b*x) " (m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~p = 0, 

where U = A*b-B*a, u = m + n + 2*p + 2, v = m + n+ 1. 



11.3. Integrands with linear-quadratic cross-degeneracy. Here the root of 
the first linear factor is required to coincide with one root of the quadratic 
whence the corresponding resultant must vanish: ra*a-rb*b = 0. To exclude ad- 
ditional degeneracies, the recurrences should be applied only if also (a*d-b*c)* 
(4*e*g-f ~2) ^0. 
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(11.3.1) u*(a*d-b*c)*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT(((A*d-B*c)*u+ (n+ 1) * (A*re - B*rf ) - V*v*c - V*v*d*x) 
*(a + b*x)~ (m+ l)*(c + d*x) ~n*(e + f *x + g*x~2)~p, x) 

- V*b* (a + b*x) ~m * (c + d*x) ~ (n + 1) * (e + f *x + g*x~2) ~ (p + 1) = 0, 

where V = A*b-B*a, u= (m + p+ l)*ra, v= (m + n + 2*p + 3)*g. 

(11.3.2) u*d*g*INT( (A+ B*x) * (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT((u*U*a+V*c + W*rf + (u*U*b + V*d + W*re)*x) 

*(a + b*x) ~ (m- l)*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
-B*b 

*(a + b*x) ~ (m- l)*(c + d*x) ~ (n+ l)*(e + f *x + g*x~2)~(p + 1) =0, 
where U= (A*d-B*c)*g, V= (m + p)*B*ra, W= (n+l)*B, u = m + n + 2*p + 2. 

(11.3.3) (n+ 1) * 1/2 * (rc*c - rd*d) *b*INT( (A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT(((n+ l)*(A*re-B*rf ) - V*(v*a- (m + p + l)*ra) - V*v*b*x) 
*(a + b*x) ~m*(c + d*x) " (n+ 1) * (e + f *x + g*x~2) ~p, x) 

- V*b*(a + b*x) ~m*(c + d*x) " (n+ 1) * (e + f *x + g*x~2) " (p+ 1) =0, 

where V = A*d - B*c, v = (m + n + 2*p + 3) *g. 

(11.3.4) (m + n + 2*p + 2)*b*g*INT((A+B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT ( (U*c + n*B*rf + (U*d + n*B*re) *x) 

*(a + b*x)~m*(c + d*x)~(n- 1) * (e + f *x + g*x~2) "p, x) 

- B*b* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) " (p + 1) = 0, 

where U= (m + n + 2*p + 2)*(A*b-B*a)*g + (m + p+ l)*B*ra. 

(11.3.5) (p+ l)*u*(rc*c-rd*d)*INT((A + B*x) 

* (a + b*x) "m * (c + d*x) "n * (e + f *x + g*x~2) ~p , x) 

- INT ( ( ( (w + n) *U*u + (n + 1) * (V - W*f ) ) *d - 2* (w + m) *W*c*g 

- 2*(w + m + n + l)*W*d*g*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+l) , x) 
+ (U*u + V - W*f - 2*W*g*x) 

*(a + b*x) ~m*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~(p+ 1) =0, 

where U = A*d-B*c, V = m* (A*se - B*sf ) , W = m*(A*re-B*rf ) - 

(m + p+ l)*(A*rc-B*rd)*b, u= (m + p + 1) * (4*e*g- f ~2) *b, w = 2*p + 3. 

(11.3.6) 2*v*(v+d)*b*g*INT((A+B*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
-INT((U*(m*rf -w*rd) + V*sf + W*w*c 
+ (U* (m*re - w*rc) + V*se + W*w*d) *x) 
*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 

- (U*b + B*(m*ra*d + v*b*f ) + 2*B*v*b*g*x) 

* (a + b*x) ~m * (c + d*x) ~ (n + 1) * (e + f *x + g*x~2) "p = 0, 

where U = 2* (m + n + 2*p + 2) * (A*d- B*c) *g + (n+l)*B*rc, V = m*(n+ l)*B*d, 
W= (n + 2*p+ l)*B*(4*e*g-f ~2)*d, v= (m + n + 2*p+ l)*d, w= (m + p)*b. 
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(11.3.7) (m + p + l)*ra*b*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
+ INT((n*(A*re-B*rf ) -V*c-V*d*x) 

*(a + b*x) ~ (m+ l)*(c + d*x) "(n- 1) * (e + f *x + g*x~2) ~p , x) 
+ (A*b-B*a)*b 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~ (p + 1) = 0, 

where V= (m + n + 2*p + 2)*(A*b-B*a)*g + (m + p+ l)*B*ra. 

(11.3.8) (n+l)*re*d*INT((A+B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- b*INT ( (U*rf + V* (v*a - (m + p) *ra) + (U*re + V*v*b) *x) 
*(a + b*x) ~ (m- l)*(c + d*x) ~ (n+ l)*(e + f *x + g*x~2)~p, x) 

+ V*b~2 

*(a + b*x) ~ (m- l)*(c + d*x) ~ (n+ l)*(e + f *x + g*x~2)~(p+ 1) =0, 
where U= (n+ 1)*B, V = A*d-B*c, v= (m + n + 2*p + 2) *g. 

(11.3.9) (p+ l)*se*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
-b*INT(((n+ l)*(A*re-B*rf) - (n + p + 2)*(A*d-B*c)*ra 

- v*V*c - v*V*d*x) 

*(a + b*x) ~ (m- l)*(c + d*x)~n*(e + f*x + g*x~2)~(p + 1) , x) 

- V* (a + b*x) ~m * (c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~ (p+ 1) = 0, 

where V = A*b*g + B*(a*g-b*f ), v = m + n + 2*p + 3. 

(11.3.10) u*b~2*d*INT((A+B*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ INT ( (u*U* (a*g - b*f ) + V*c - W*rf - (u*U*b*g - V*d + W*re) *x) 
*(a + b*x)~ (m+ l)*(c + d*x)~n*(e + f *x + g*x~2) ~ (p- 1) , x) 

- B*b* (a + b*x) ~ (m + 1) * (c + d*x) " (n + 1) * (e + f *x + g*x~2) ~p = 0, 

where U = A*d-B*c, V= (n + p+ l)*B*ra, W= (n+ 1)*B, u = m + n + 2*p + 2. 

(11.3.11) 2*(m + p + l)*(p + l)*(4*e*g-f ~2) *b*INT( (A + B*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT((n*U*d+ v*V*c+ (v + n)*V*d*x) 

*(a + b*x)~m*(c + d*x) ~(n- l)*(e + f *x + g*x~2) ~ (p+ 1) , x) 
+ (U + V*x)*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+ 1) =0, 

where U = m*A*ra+ (m + 2*p + 2) * (A*f - 2*B*e) *b, 

V= (m+2*p + 2)*(2*A*g-B*f )*b + m*B*ra, v = m + 2*p + 3. 

(11.3.12) (m + n + 2*p + 2)*v*b*d*INT((A+B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
+ INT((U*(u*f -w) -B*u*v*e + (U*u*g-B*v*w)*x) 

*(a + b*x)~m*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~ (p- 1) , x) 
-b*(U + B*v*x) 

*(a + b*x)~m*(c + d*x)~(n + 1) * (e + f *x + g*x~2) ~p = 0, 

where U= (m + n + 2*p + 2)*A*d- (m + 2*p+ l)*B*c, u= (m + 2*p)*b, v= (n+ l)*d, 
w = m*a*g + p*b*f . 
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11.4. Integrands with a self-degenerate quadratic. Here the quadratic poly- 
nomial is required to possess a double root whence its discriminant must vanish: 
4*e*g-f~2 = 0. To exclude additional degeneracies, the recurrences should be ap- 
plied only if also ra*rc + 0. 

(11.4.1) u*(a*d-b*c)*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT(((A*d-B*c)*u- V*(v*c- (n+ l)*rc) - V*v*d*x) 
*(a + b*x) " (m+ l)*(c + d*x) ~n * (e + f *x + g*x~2) ~p, x) 

-V*(f + 2*g*x) 

* (a + b*x) " (m+ 1) * (c + d*x) " (n+ 1) * (e + f *x + g*x~2) ~p = 0, 
where V = A*b-B*a, u= (m+ l)*ra, v = 2*(m + n + 2*p + 3)*g. 

(11.4.2) 2* (m + n + 2*p + 2) *d*g*INT ( (A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
-INT((U*a + V*f + (U*b + 2*V*g)*x) 

*(a + b*x) ~ (m- l)*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
-B*(f + 2*g*x)*(a + b*x) ~m*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~p = 0, 

where U = 2* (m + n + 2*p + 2) * (A*d- B*c) *g + (m + n+ l)*B*rc, 
V = m*B*(a*d-b*c). 

(11.4.3) u*w*INT ( (A + B*x) * (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 
+ INT( (v* (U*a*d + V*b*c - A*w*b*d) 

- (U + V-B*w)*((v + m+ l)*b*c + (v + n+ l)*a*d) 

- (v + m + n + 2)*(U + V-B*w)*b*d*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f *x + g*x~2) ~ (p+ 1) , x) 
-u*(A*f -2*B*e + (2*A*g-B*f )*x) 

* (a + b*x) ~ (m + 1 ) * (c + d*x) ~ (n + 1 ) * (e + f *x + g*x~2) ~p 
+ (U + V-B*w) 

* (a + b*x) ~ (m + 1) * (c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) " (p + 1) = 0, 

where U= (m + 2*p + 2) * (A*rc - B*rd) *b, V= (n + 2*p + 2)*(A*ra-B*rb)*d, 
u= (p + l)*(ra*c-rb*d), v = 2*(p+ 1), w = (2*p + 1) * (ra*c - rb*d) . 

(11.4.4) 2*t*(t + b*d)*INT((A + B*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ INT ( (V*u + v* (V*b + W*a) *rd + w* (V*d + W*c) *rb 
+ (v* (V*b + W*a) *rc + w* (V*d + W*c) *ra - W*u) *x) 
*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 

- (V*f + 2*W*e + (2*V*g + W*f )*x) 

*(a + b*x)~(m+ l)*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~ (p - 1) 
-2*B*t*(a + b*x)~(m+ l)*(c + d*x)~(n + l)*(e + f *x + g*x~2)~p = 0, 

where V= (m + n + 2*p + 2)*A*b*d-B*((m+2*p + l)*b*c + (n + 2*p+ l)*a*d), 
W = 2*p*B*b*d, t= (m + n + 2*p + l)*b*d, u= (2*p- 1) * (ra*c - rb*d) , v = m + 2*p, 
w = n + 2*p. 

(11.4.5) u*b*INT((A + B*x)*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- INT(((B*u + V*v)*c-n*V*rc+ (B*u + V*v)*d*x) 

*(a + b*x)~(m+ l)*(c + d*x) ~(n- l)*(e + f *x + g*x~2)~p, x) 
+ V* (f + 2*g*x) * (a + b*x) ~(m+l)*(c + d*x) ~n * (e + f *x + g*x~2) "p = 0, 

where V = A*b-B*a, u= (m+ l)*ra, v = 2*(m + n + 2*p + 2)*g. 
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(11.4.6) (2*p+ l)*u*rc*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- INT((V*(v*a*d-m*w) + W*(v*b*c + (n+l)*w) + v*(V + W)*b*d*x) 
*(a + b*x)~(m- l)*(c + d*x)~n*(e + f*x + g*x~2)~(p + 1) , x) 

-u*(A*f - 2*B*e+ (2*A*g-B*f )*x) 

* (a + b*x) ~m * (c + d*x) ~ (n + 1 ) * (e + f *x + g*x~2) "p 

+ (V + W)*(a + b*x)~m*(c + d*x)~(n+ l)*(e + f*x + g*x~2)~(p + 1) =0, 

where V= (n + 1) * (A*ra- B*rb) *d + (2*p + l)*(A*d-B*c)*ra, 

W = m*(A*rc-B*rd)*b, u= (p+ 1) * (ra*c - rb*d) , v = m + n + 2*p + 3, w = a*d-b*c. 

(11.4.7) (m + n + 2*p + 2) *v*w*b*d*INT ( (A + B*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ INT((U*((n+ l)*rb*d+ (2*p- l)*ra*c) + (U + p*B*v)*w*rd 
+ (U*(n + 2*p)*ra*d + (U + p*B*v)*w*rc)*x) 
*(a + b*x)~(m+l)*(c + d*x)~n*(e + f*x + g*x~2)~(p- 1) , x) 
-U*(rb + ra*x) 

*(a + b*x)~ (m+ l)*(c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) " (p - 1) 
-B*v*w*(a + b*x)~(m+ 1) * (c + d*x) " (n+ 1) * (e + f *x + g*x~2) ~p = 0, 

where U= (m+ l)*(A*d-B*c)*b + (n + 2*p + 1) * (A*b - B*a) *d, v = 2*(a*d-b*c), 
w= (m+ l)*b. 

11.5. Doubly degenerate integrands. Now the roots of both linear factors and 
one root of the quadratic are required to coincide whence the following must van- 
ish: a*d - b*c = ra*a - rb*b = 0. To exclude additional degeneracies, the recurrences 
should be applied only if also ra^ 0. Note that each relation can be used backwards 
for its own inverse. 

(11.5.1) u*ra*INT((a+b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- (u + p+ l)*g 

*INT((a + b*x)~(m+ 1) * (c + d*x) ~n* (e + f *x + g*x~2) "p , x) 
+ b* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~ (p+ 1) = 0, 

where u = m + n + p + l. 

(11.5.2) (m + n + p + l)*(p+ l)*(4*e*g-f ~2)*b 

*INT((a + b*x)~m*(c + d*x)~n*(e + f *x + g*x~2)~p, x) 

- w*(w + l)*b*g 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+l) , x) 

- ( (p + 1) *ra- w*a*g - w*b*g*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+ 1) =0, 

where w = m + n + 2*p + 2. 

(11.5.3) b*INT((a + b*x)~m*(c + d*x)~n*(e + f *x + g*x~2)~p, x) 

- d*INT((a+b*x)~ (m+ l)*(c + d*x)~(n- l)*(e + f *x + g*x~2)~p, x) 
= 0. 

(11.5.4) (p + l)*ra*INT((a+b*x) ~m*(c + d*x) ~n*(e + f *x + g*x~2)~p, x) 
+ (m + n + 2*p + 2)*b~2 

*INT((a + b*x)~(m- 1) * (c + d*x) ~n* (e + f *x + g*x~2) " (p+ 1) , x) 

- b* (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~ (p + 1) = 0. 

Now the roots of the two linear factors are required to coincide while the qua- 
dratic is required to possess a double root whence the following must vanish: 
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a*d - b*c = 4*e*g - f ~2 = 0. To exclude additional degeneracies, the recurrences 
should be applied only if also ra^O. Note that each relation can be used back- 
wards for its own inverse. 

(11.5.5) (m + n+ l)*ra*INT((a+b*x) ~m*(c + d*x) ~n * (e + f *x + g*x~2) ~p, x) 
-2*(m + n + 2*p + 2)*g 

*INT((a + b*x)~(m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~p, x) 
+ (f + 2*g*x)*(a + b*x)~(m+ l)*(c + d*x) ~n*(e + f *x + g*x~2) ~p = 0. 

(11.5.6) (p+ l)*(2*p + l)*(ra*a-rb*b) 

*INT( (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- w*(w + b) 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+l) , x) 

- (p+ l)*(rb + ra*x) 

* (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~p 
+ w*(a + b*x)~(m + 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~ (p+ 1) =0, 

where w = (m + n + 2*p + 2) *b. 

(11.5.7) b*INT((a + b*x)~m*(c + d*x)~n*(e + f *x + g*x~2)~p, x) 
-d*INT((a + b*x)~(m + 1) * (c + d*x) ~ (n- 1) * (e + f *x + g*x~2) ~p, x) 
= 0. 

(11.5.8) (p+l)*(2*p+l)*ra 

*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- (m + n + 2*p + 2)*w*b 

*INT((a+b*x)~(m- 1) * (c + d*x) ~n * (e + f *x + g*x~2) " (p+ 1) , x) 

- (p + 1) * (rb + ra*x) * (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p 
+ w*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p + 1) =0, 

where w = (m + n) *b. 

Now the roots of the two linear factors are required to coincide with different 
roots of the quadratic whence the following must vanish: re = rf = 0. To exclude 
additional degeneracies, the recurrences should be applied only if also a*d - b*c ^ 0. 
Note that each relation can be used backwards for its own inverse. 

(11.5.9) (m + p+ l)*(a*d-b*c) 

*INT( (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) ~p , x) 

- (m + n + 2*p + 2)*d 

*INT((a + b*x)~(m+ 1) * (c + d*x) ~n * (e + f *x + g*x~2) ~p, x) 
+ (a + b*x) ~ (m+ 1) * (c + d*x) ~ (n+ 1) * (e + f *x + g*x~2) ~p = 0. 

(11.5.10) u*v*(ra*c-rb*d) 

*INT((a+b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- (u + v)*(u + v + l)*b*d 

*INT((a+b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p + 1) , x) 
+ (u*a*d+ v*b*c + (u + v)*b*d*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p + 1) =0, 

where u = m + p+ l, v = n + p + 1 . 

(11.5.11) (m + p+ l)*b*INT((a+b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ (n + p)*d 

*INT((a+b*x) " (m+ 1) * (c + d*x) " (n- l)*(e + f *x + g*x~2) ~p, x) 

- (a + b*x)~(m+ l)*(c + d*x)~n*(e + f*x + g*x~2)~p = 0. 
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(11.5.12) (n + p+ l)*ra*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ (m + n + 2*p + 2)*b~2 

*INT((a + b*x)~(m- l)*(c + d*x) ~n*(e + f *x + g*x~2) " (p+ 1) , x) 
-b*(a + b*x)~m*(c + d*x)~n*(e + f *x + g*x~2) ~ (p + 1) =0. 

Now the quadratic is required to possess a double root which coincides with 
the root of the first linear factor whence the following must vanish: ra = rb = 0. 
To exclude additional degeneracies, the recurrences should be applied only if also 
a*d-b*c^0. Note that each relation can be used backwards for its own inverse. 

(11.5.13) u*(a*d-b*c)*INT((a+b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- (u + n+ l)*d 

*INT((a + b*x)~ (m+ l)*(c + d*x) ~n*(e + f *x + g*x~2) ~p, x) 
+ (a + b*x) ~ (m + 1) * (c + d*x) " (n + 1) * (e + f *x + g*x~2) ~p = 0, 

where u = m + 2*p + 1. 

(11.5.14) (n+ l)*(a*d-b*c) 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ (m + n + 2*p + 2)*b 

*INT((a + b*x)~m*(c + d*x)~(n + l)*(e + f *x + g*x~2) ~p, x) 

- (a + b*x) ~ (m + 1) * (c + d*x) " (n + 1) * (e + f *x + g*x~2) ~p = 0. 

(11.5.15) (m + 2*p + l)*u*l/2*(rc*c-rd*d) 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- v*(v + d) 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~(p+l) , x) 
+ (u*(a*d-b*c) + v*a + v*b*x) 

*(a + b*x)~ (m- l)*(c + d*x)~ (n+ l)*(e + f *x + g*x~2)~ (p+ 1) =0, 

where u = m + 2*p + 2, v=(m + n + 2*p + 2) *d. 

(11.5.16) (m + 2*p + l)*b*INT((a + b*x)~m ! Kc + d*x)~n*(e + f*x + g*x~2)~p, x) 
+ n*d 

*INT((a + b*x)~(m+ l)*(c + d*x) ~(n- l)*(e + f *x + g*x~2)~p, x) 

- (a + b*x)~(m + l)*(c + d*x)~n*(e + f *x + g*x~2) ~p = 0. 

(11.5.17) u*rc*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x + g*x~2) "p , x) 
+ 2*(u + n+ l)*b*d 

*INT((a+b*x)~(m- l)*(c + d*x) ~n*(e + f *x + g*x~2) " (p+ 1) , x) 
-2*b*(a + b*x)~ (m- 1) * (c + d*x) " (n+ 1) * (e + f *x + g*x~2) " (p+ 1) 
= 0, 

where u = m + 2*p + 1. 

(11.5.18) (m + 2*p+ l)*u*rc 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x + g*x~2)~p, x) 

- 2*n*v*d 

*INT((a + b*x)~m*(c + d*x)~ (n- 1) * (e + f *x + g*x~2) ~ (p+ 1) , x) 

- 2* (u*b*c - v*a - n*b*d*x) 

*(a + b*x)~ (m- l)*(c + d*x)~n*(e + f *x + g*x~2) ~ (p+ 1) =0, 



where u = m+2*p + 2, v= (m + n + 2*p + 2)*d. 
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12. Integrands involving (a + b*x)~m*(c + d*x)~n*(e + f*x)~p*(g + h*x)~q 
Abbreviations used with these integrands: 

ra= (e*h + f *g) *b - a*f *h, rb = b*e*g, re = (a*h + b*g) *f -b*e*h, rf = a*f *g. 
saa = b*c*f*h + ra*d, sab = a*c*f*h + ra*c + rb*d, sbb = rb*c. 

12.1. Nondegenerate integrands. To exclude integrands with confluent roots, 
the following recurrences should be applied only if the overall discriminant 
does not vanish: (a*d - b*c) * (a*f - b*e) * (a*h- b*g) * (c*f - d*e) * (c*h- d*g) * 
(e*h-f*g) *0. 

(12.1.1) (m+ l)*(a*d-b*c)*(a*f -b*e)*(a*h-b*g)*INT((A + B*x + C*x~2) 

* (a + b*x) ~m * (c + d*x) ~n* (e + f *x) ~p * (g + h*x) ~q, x) 
+ INT(((m+l)*(A*(saa*a-sab*b) + (B*b - C*a) *sbb) 

-W*((n+ l)*d*e*g+ (p+ l)*c*f *g+ (q+ l)*c*e*h) 

- ((m+ l)*((A*b-B*a)*saa + C*(sab*a-sbb*b)) 

+ W*((n + p + 2)*d*f*g + (n + q+2)*d*e*h + (p + q+ 2) *c*f *h) ) *x 

- (m + n + p + q + 4)*W*d*f*h*x~2) 

*(a + b*x) ~ (m+ l)*(c + d*x) ~n*(e + f*x) ~p*(g + h*x) ~q, x) 
+ W*(a + b*x) ~(m + l)*(c + d*x) ~ (n+ 1) * (e + f *x) ~ (p+ 1) 
*(g + h*x)- (q+1) =0, 

where W = A*b~2-B*a*b + C*a~2. 

(12.1.2) u*INT((A + B*x + C*x~2)*(a+b*x)-m 

* (c + d*x) ~n * (e + f*x) ~p * (g + h*x) ~q, x) 

- INT(((A*u-V)*a-m*C*sbb + ((A*u-V)*b 

+ (B*u-W)*a-m*C*sab)*x + ( (B*u - W) *b - m*C*saa) *x~2) 
*(a + b*x) ~ (m- l)*(c + d*x)~n*(e + f*x)~p*(g + h*x)~q, x) 

-C*(a + b*x) ~m*(c + d*x) ~(n + l)*(e + f *x) " (p+ l)*(g + h*x) " (q+ 1) 

= 0, 

where V = C*((n + l)*d*e*g+ (p+ l)*c*f*g + (q+ l)*c*e*h), 
W = C*((n + p + 2)*d*f *g+ (n + q + 2)*d*e*h + (p + q+ 2) *c*f *h) , 
u= (m + n + p + q + 3)*d*f *h. 

(12.1.3) u*b*INT((A + B*x + C*x'2) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

+ INT((U*c + V*(n*d*e*g + (p+ l)*c*f *g + (q+ l)*c*e*h) + (U*d 
-C*u*c + V*(v*c*f*h + (n + p+ l)*d*f *g+ (n + q+ l)*d*e*h))*x 

- (C*u- (v + n)*V*f *h)*d*x~2) 

*(a + b*x) " (m+ l)*(c + d*x)~(n- l)*(e + f *x) ~p*(g + h*x) ~q, x) 

- V*(a + b*x) ~(m+ 1) * (c + d*x) "n * (e + f *x) ~ (p+ 1) * (g + h*x) " (q+ 1) 
= 0, 

where U= (m+ l)*(A*ra*b- (B*b - C*a) *rb) , V = A*b~2- (B*b-C*a)*a, 
u= (m+l)*(a*f -b*e)*(a*h-b*g), v = m + p + q + 3. 

12.2. Singly degenerate integrands. Here the roots of the first two linear fac- 
tors are required to coincide whence the corresponding resultant must vanish: 
a*d-b*c = 0. To exclude additional degeneracies, the recurrences should be ap- 
plied only if also (a*f - b*e) * (a*h- b*g) * (e*h- f *g) + 0. 
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(12.2.1) u* (a*f - b*e) * (a*h - b*g) *INT( (A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f*x) ~p * (g + h*x) ~q, x) 
+ INT((u*(A*ra-B*rb) + V*((q+ l)*e*ri + (p+l)*f*g) 

+ (u + p + q + 2)*V*f*h*x) 

* (a + b*x) " (m+ 1) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 
-V*(a + b*x)~(m+l)*(c + d*x)~n*(e + f*x)~(p+l)*(g + h*x)~(q+l) 
= 0, 

where V = A*b-B*a, u = m + n+ 1. 

(12.2.2) (m + n + p + q + 2)*f*h*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

- INT ( (U*a - V*rb + (U*b - V*ra) *x) 

*(a + b*x) ~ (m- l)*(c + d*x)~n*(e + f*x)~p*(g + h*x)~q, x) 

- B* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) " (p + 1 ) * (g + h*x) ~ (q+ 1) = 0, 

where U= (m + n + p + q + 2)*A*f *h-B*((p + l)*f *g + (q+ l)*e*h), V= (m + n)*B. 

(12.2.3) (p + l)*(a*f -b*e)*(e*ri-f*g)*INT((A + B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

- INT(((p+ l)*(A*re-B*rf ) + V*((m + n+ l)*b*g + (q+ l)*a*ri) 

+ (m + n + p + q + 3) *V*b*h*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~ (p + 1) * (g + h*x) "q, x) 

+ V*(a + b*x)~(m+l)*(c + d*x)"n*(e + f*x)~(p+l)*(g + h*x)"(q+l) 
= 0, 

where V = A*f - B*e. 

(12.2.4) (m + n + p + q + 2)*b*h*INT((A + B*x) 

* (a + b*x) ~m* (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

- INT ( (U*e - p*B*rf + (U*f - p*B*re) *x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x)~(p- l)*(g + h*x)~q, x) 

- B* (a + b*x) ~ (m + 1) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) " (q+ 1) = 0, 

where U = A*(m + n + p + q+2)*b*h-B*((m + n+ l)*b*g+ (q+ l)*a*ri). 

(12.2.5) b*INT((A+B*x) 

* (a + b*x) ~m * (c + d*x) "n * (e + f *x) ~p * (g + h*x) ~q, x) 

- d*INT((A + B*x)*(a+b*x)~ (m + 1) * (c + d*x) " (n- l)*( e + f *x)~p 
*(g + h*x)~q, x) =0. 

(12.2.6) u*b*INT((A+B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

- INT((B*u*e + V*(v*e*h + p*f *g) + (B*u + (v + p)*V*h)*f *x) 

*(a + b*x)~(m+ l)*(c + d*x) ~n*(e + f *x)~(p- l)*(g + h*x) ~q, x) 
+ V* (a + b*x) ~ (m + 1) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~ (q+ 1) = 0, 

where V = A*b-B*a, u= (m + n+ l)*(a*h-b*g), v = m + n + q+2. 

(12.2.7) u*f*INT((A+B*x) 

*(a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

- INT( (B*u*a + V* (v*a*h + (m + n)*b*g) + (B*u + (v + m + n)*V*h)*b*x) 
*(a + b*x)~(m- l)*(c + d*x) ~n*(e + f *x)~(p + l)*(g + ri*x) ~q, x) 

+ V* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~ (p + 1 ) * (g + h*x) ~ (q+ 1) = 0, 

where V = A*f -B*e, u= (p+ l)*(e*h-f*g), v = p + q + 2. 
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(12.2.8) u*f*INT((A+B*x) 

* (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 
-INT((B*u*g-V*(v*b*g+q*a*h) + (B*u- (v + q) *V*b) *h*x) 

*(a + b*x)~m*(c + d*x)~n*(e + f*x)~(p + l)*(g + h*x)~(q- 1) , x) 

- V* (a + b*x) ~ (m+ 1) * (c + d*x) "n * (e + f *x) ~ (p + 1) * (g + h*x) ~q = 0, 

where V = A*f -B*e, u= (p + 1) * (a*f - b*e) , v = m + n + p + 2. 

12.3. Doubly degenerate integrands. Now the roots of the first three linear 
factors are required to coincide whence the following must vanish: a*d-b*c = 
a*f-b*e = 0. To exclude additional degeneracies, the recurrences should be ap- 
plied only if also a*h-b*g^ 0. Note that each relation can be used backwards for 
its own inverse. 

(12.3.1) u*(a*h-b*g) 

*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) "p * (g + h*x) ~q , x) 

- (u + q + l)*h 

*INT( (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) "q, x) 
+ (a + b*x) ~ (m+ l)*(c + d*x)~n*(e + f *x) ~p*(g + h*x) ~(q+ 1) =0, 

m + n + p + 1. 

(q+ l)*(a*h-b*g) 

*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q , x) 
+ (m + n + p + q+2)*b 

*INT((a + b*x)~m*(c + d*x)~n*(e + f*x)~p*(g + h*x)~(q + 1) , x) 

- (a + b*x) ~ (m+ l)*(c + d*x)~n*(e + f *x) ~p*(g + h*x) ~ (q+ 1) =0. 

b*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 
-d*INT((a + b*x)~(m + 1) * (c + d*x) ~ (n- l)*(e + f*x)~p 
*(g + h*x)-q, x) =0. 

(m + n + p+ l)*b* 

INT( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 
+ q*h*INT ( (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x) ~p 
*(g + h*x)- (q- 1) , x) 

- (a + b*x) ~ (m+ 1) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q= 0. 

Now the roots of the first two linear factors are required to coincide as are the 
roots of the last two linear factors whence the following must vanish: a*d-b*c = 
e*h - f *g = 0. To exclude additional degeneracies, the recurrences should be applied 
only if also a*f - b*e # 0. Note that each relation can be used backwards for its own 
inverse. 

(12.3.5) (m + n+l)*(a*f -b*e) 

*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q , x) 

- (m + n + p + q+2)*f 

*INT((a + b*x)~(m+ 1) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 
+ (a + b*x)~(m+l)*(c + d*x)~n*(e + f*x)~(p+l)*(g + n*x)~q=0. 

(12.3.6) b*INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 

- d*INT((a + b*x)~ (m+ 1) * (c + d*x) " (n- l)*(e + f *x)~p 
*(g + h*x)~q, x) =0. 
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(12.3.7) (m + n+l)*b* 

INT ( (a + b*x) ~m * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q, x) 
+ (p + q)*f *INT((a+b*x)~ (m+ l)*(c + d*x) ~n*(e + f *x)~(p- 1) 

*(g + h*x) ~q, x) 
- (a + b*x) " (m + 1 ) * (c + d*x) ~n * (e + f *x) ~p * (g + h*x) ~q = 0. 
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